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We investigate the current-induced spin polarization in the two-dimensional hole gas (2DHG) with 
the structure inversion asymmetry. By using the perturbation theory, we re-derive the effective k- 
cubic Rashba Hamiltonian for 2DHG and the generalized spin operators accordingly. Then based 
on the linear response theory we calculate the current-induced spin polarization both analytically 
and numerically with the disorder effect considered. We have found that, quite different from the 
two-dimensional electron gas, the spin polarization in 2DHG depends linearly on Fermi energy in 
the low doping regime, and with increasing Fermi energy, the spin polarization may be suppressed 
and even changes its sign. We predict a pronounced peak of the spin polarization in 2DHG once the 
Fermi level is somewhere between minimum points of two spin-split branches of the lowest light-hole 
subband. We discuss the possibility of measurements in experiments as regards the temperature 
and the width of quantum wells. 

PACS numbers: 72.25.-b, 85.75.-d, 71.70.Ej, 72.25.Pn 



I. INTRODUCTION 

In order to reduce the electric leakage and to meet the 
challenge brought about by the reduced physical size of 
the future nano-electronics, it is being explored to re- 
place the electron charge with the spin degree of freedom 
in the electronic transport. This is the ambitious goal 
of researchers in the field of spintronics i 1 * 2 * 3 - One of basic 
issues in this field is how to generate the polarized spin 
in devices. As an straightforward way, the spin injection 
from ferromagnetic layers may provide a possible solution 
to this problem if the interface mismatch problem can be 
avoided, but it is more desirable to generate spin polar- 
ization directly by electric means in devices because of 
its easy controllability and compatibility with the stan- 
dard microelectronics technology. 1 * 2 ^ The spin-orbit cou- 
pling (SOC) in semiconductors, which relates the elec- 
tron spin to its momentum, may provide a controllable 
way to realize such purpose. Based on this idea, the phe- 
nomenon of current-induced spin polarization (CISP) has 
recently attracted extensive attentions of a lot of research 
groups 4- 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 ' 14 - 15 ' 16 - 17 ' 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 

As early as in 1970's, the CISP due to the spin-orbit 
scattering near the surface of semiconductor thin films 
was predicted by Dyakonov and Pereli^ Restricted by 
experimental conditions at that time, this prediction 
was ignored until the beginning of 1990's. With the 
development of sample fabrication and characterization 
technology in low-dimensional semiconductor systems, it 
was realized that such phenomena could also exist in 
quantum wells and heterostructures with the structure 
or bulk inversion asymmetry£>£ Later, many interest- 
ing topics about CISP have been raised, such as the 
joint effect of the Rashba and Dresselhaus SOC mech- 
anism,— vertex correctional 7 ^ quantum correction 9 ^ 
and resonant spin polarization^ Experimentally, CISP 
was first observed by Silov et al 12 in two-dimensional 



hole gas (2DHG) by using the polarized photolumines- 
cence] 27 : 28 ' 29 When inputting an in-plane current into 
the 2DHG system, they observed a large optical polar- 
ization in photoluminescence spectral Later, Kato et 
al demonstrated the existence of the CISP in strained 
nonmagnetic semiconductors ; 13 i 14 and Sih et al detected 
the CISP in the two-dimensional electron gas (2DEG) 
in (110) AlGaAs quantum well^ The CISP was also 
found in ZnSe epilayers even up to the room temper- 
ature^ Very recently, the converse effect of CISP has 
been clearly shown by Yang et al experimentallyji 7 - and 
the spin photocurrent has also been observed. 18 : 19 ' 20 

So far most theoretic investigations about 
the CISP deal with the electron SOC sys- 
tems, 4 ! 5 ! 6 ! 7 ! 8 ! 9 . 10 ! 21 ! 22 ! 23 : 24 ! 25 . 26 Thus the CISP in 
the 2DHG system as shown in Silov's experiments 
was also interpreted in terms of the linear-fc Rashba 
coupling of the 2DEG systems with several parameters 
adjusted*^ As we shall show later, this treatment is 
not appropriate for 2DHG. Unlike the electron system, 
the hole state in the Luttinger-Kohn Hamiltonian^ is 
a spinor of four components. As each component is 
a combination of spin and orbit momentum, the spin 
of a hole spinor is not a conserved physical quantity. 
Therefore, the " spintronics" for hole gas is in fact a 
combination of spintronics and orbitronics 3 - 1 -. If only 
the lowest heavy hole (HH1) subband is concerned, by 
projecting the multi-band Hamiltonian of 2DHG with 
structural inversion asymmetry into a subspace spanned 
by | ± |) mostly relevant with the HH1 states, we 
can obtain the fc-cubic Rashba mode l 32 : 33 ' 34 : 35 ' 3 ^. We 
emphasize here in this lowest heavy hole subspace, the 
spin operators are no longer represented by three Pauli 
matrices, because the "generalized spin" we shall adopt 
is a hybridization of spin and orbit angular momentum. 
In deriving the effective Hamiltonian from the Luttinger- 
Kohn Hamiltonian by the perturbation and truncation 
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procedure to higher orders, one must take care of the 
corresponding transformation for the spin operator in 
order to obtain the correct expression. In the following, 
we will use the terminology " generalized spin" , or the 
"spin" for short, to denote the total angular momentum 
in the spin-orbit coupled systems. 

The aim of the present paper is to investigate the CISP 
of 2DHG in a more rigorous way. Namely, we will derive 
the A;-cubic Rashba model and the corresponding spin op- 
erators for holes, and on this basis we will present both 
analytical and numerical results for the CISP in 2DHG. 
This paper is organized as follows. In Sec II the gen- 
eral formalism and the Hamiltonian for the 2DHG with 
structural inversion asymmetry is given. In Sec III in 
the low doping regime, with the perturbation theory, the 
Hamiltonian and spin operators in the lowest heavy hole 
subspace are derived, and applied to analytical calcula- 
tion of the CISP in 2DHG. In Sec IV, we will show the 
numerical calculations agree well with the analytical re- 
sults at the low-doping regime; while in the high doping 
regime the numerical results predict some new features of 
CISP. Particularly, we predict a pronounced CISP peak 
when Fermi energy lies little above the energy minimum 
of the lowest light hole (LH1) subband. Finally, a brief 
summary is drawn. 
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where 71,72 is the Luttinger parameters, mo is the free 
electron mass, the in-plane wave vector k = (k x , k y ), de- 
noted in the polar coordinate as k = (k, 9), \i± = k x ±ik y 
and k z = —id/dz. The other terms, such as anisotropic 
term, C terms or hole Rashba term. ) 35 i 36 i 37 have only neg- 
ligible effects and are omitted in our calculation. Cor- 
respondingly, the x-, y-, z- component of the "spin"-! 
operator respectively reads 
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II. FORMALISM 



A. Hole Hamlitonian 



A p-doped quantum well system with structural in- 
version asymmetry can be described as the isotropic 
Luttinger-Kohn Hamiltonian with a confining asymmet- 
rical potential, 



H = H L + V c (z) + V a (z). 



(1) 



Here in order to compare the analytical results with 
the numerical one, the confining potential along the z- 
direction V c (z) is taken as 



Vc(z) 



-L z /2 < z < L z /2 
00 otherwise, 



(2) 



where L z is the well width of the quantum well. The 
asymmetrical potential, which stems from a build-in elec- 
tric field F via the gate voltage or <5-doping is V a {z) — 
eFz, which breaks the inversion symmetry and lifts the 
spin doublet degeneracy. 

Let S be the generalized spin operator of a hole state, 
and S z be the z-component of S, the isotropic Luttinger- 
Kohn Hamiltonian Hl in the |5, S z ) representation (four 
basis kets written in the sequence of {||), | — | — 
|)}) is expressed as 
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We stress here again that the " spin" of the | spinor is ac- 
tually its total angular momentum, which is a linear com- 
bination of spin and orbit angular momentum of a valence 
band electron. In polarized optical experiments, such as 
polarized photoluminescenc o 27 i 28 ' 29 or Kerr/Farady rota- 
tio n 13 i 14 , it is appropriate to introduce such a generalized 
spin. 

For the infinitely confining potential, we expand the 
eigenfunction 0„ associated with the vth hole subband 
in terms of confined standing waves as 
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with 



> > /2 fnir(z + L z /2)\ . 

|n ' h) = yT z s ( I z ) 1 ( } 

where r = (x,y), n is the confinement quantum 
number for the standing wave along the z-direction, 
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and A/, denotes the A^-component of the hole (A/, = 
3/2,1/2,-1/2,-3/2). Since we are only interested in 
the low energy physics, a finite number of n will result 
in a reasonable accuracy and the effective Hamiltonian 
is reduced into a square matrix with a dimension of An. 
In this way we obtain the hole subband structure analyt- 
ically or numerically. 



B. Expression for CISP 

In the framework of the linear response theory, the 
electric response of spin polarization in a weak external 
electric field E can be formulated as 11 



(S a ) 



13 



(13) 



where (S a ) is the thermodynamically averaged value of 
the spin density. The electric spin susceptibility Xa/3 can 
be calculated by Kubo formula.— By the Green function 
formalism, the Bastin version of Kubo formula 39 reads 



dG h 
~dE 



dG 
dE 



Xati = S/ dE f( E ) Tr ( s ° 

(14) 

where G R and G A are the retarded and advanced Green 
function, respectively, A = i(G R — G A ) is the spectral 
function, f(E) is the Fermi distribution function, vp 
is the velocity operator along the (3 direction, and the 
bracket (• • -) c represents the average over the impurity 
configuration. 

To taken the vertex correction into account, we use the 
Streda-Smrcka division of Kubo formula , 39 ' 40 

Xa/3 = _ £ / dE^^-Tr{S a G R {E F )v p G A {E F )) c , 

in which we retain only the non-analytical part, and ne- 
glect the analytical part, because the latter is much less 
important in the present case. In the following, we will 
use Eq. (fl~5|) to analytically calculate the electric spin sus- 
ceptibility (ESS) with the vertex correction considered; 
meanwhile we will carry out the numerical calculation 
with Eq. p4p in the relaxation time approximation. We 
shall show that the analytical and numerical results are 
in good agreements with each other in the regime of low 
hole density. 



([T]) is invariant under the space inversion transformation 

(16) 



x — > —x, y — > — y, z — > —z, 
S x > S x , Sy > S y , S z ► S z . 



if the origin point of z-axis is set at the mid-plane of the 
quantum well. Applying the space inversion transforma- 
tion (JTljJ to Eq. (JT3JI, we have 

(Sot) = XafjEp — > (S a ) — ~XapE(i, (17) 

whereby Xap = —Xu(3- This implies that no CISP ap- 
pears when the inversion symmetry exists in the system. 
So the asymmetrical potential V a is crucial for the CISP. 

In the presence of an asymmetrical potential V a , the 
Hamiltonian JTJ) is invariant versus the rotation along z- 
axis with ^ in both the real space and the spin space, 



x -> y, y^> -x, z -> z, 

S X * Sy, Sy *■ S X , S Z * S Z 



(18) 



With the above transformations (TTg|) . Eq. I|13p will give 

(Sx) = XxyE y -> (S y ) = ~XxyE x , (19) 

(S X ) — XxxE x — > (Sy) — XxxEy. (20) 

Combined with (S y ) = XyxE x and (S y ) = XyyE y , we get 



Xxy Xyx , 

Xxx — Xyyi 



(21) 
(22) 



which are direct consequence of the rotation symmetry 
along the z-axis. 



III. ANALYTICAL RESULTS FOR CISP IN 
2DHG 

In the low hole density regime an effective Hamil- 
tonian can be obtained by projecting the Hamiltonian 
(PT|> into the subspace spanned by the lowest heavy hole 
states, which, by using the truncation approximation and 
projection perturbation method ) 35 ' 36 ' 37 ' 41 ' 42 ' 43 ' 44 ' 45 is re- 
duced to the widely used fc-cubic Rashba model. More 
importantly, the corresponding spin operators in the sub- 
space will be obtained properly, and the ESS of 2DHG 
with the impurity vertex correction will be worked out. 
Then we will compare and contrast the different behav- 
iors of the CISP in the 2DEG and 2DHG in this Section. 



C. Symmetry 



fc-cubic Rashba Model 



The general properties of Xa/3 will be critically de- 
termined by symmetry of the system. For the two- 
dimensional system we investigate, the index a(/3) in 
Eq. (|13[) is simply chosen to be x or y in the following. 
Without the asymmetrical potential V a , the Hamiltonian 



To obtain an approximate analytical expression, we 
take the following procedure. First we expand a hole 
state in terms of 8 basis wave functions associated with 
\n,X h ) (n = 1, 2 and A fc = §,£,-£,-§) ( Eq.H^. Then 
for a given k, we may express the Hamiltonian (fT|) in 
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terms of an 8 x 8 matrix, which by the perturbation proce- 
dure can be further projected into the subspace spanned 
by the |1, §) and |1, — |) states. Thus we obtain a 2 x 2 
matrix as ( See Appendix A for details), 

h 2 k 2 



H k 3 = 



2rrih 



ia(k_a+ — fci_<7_), 



(23) 



where the Pauli matrix a± = \{p x ± ia v ), the effective 



mass is renormalized into 



m h = m 7i + 72 



256 72 2 



37r 2 (3 7 i + 107 2 ) 
and the fc-cubic Rashba coefficient 



9?r 6 (37i + 107 2 )(7i - 2 72 



(24) 



(25) 



Note that Eq. (|23|) is just the fc-cubic Rashba model, in 
which the Rashba coefficient a is proportional to asym- 
metrical potential strength F, in agreement with the re- 
sults by Winkler^ We can rewrite the fc-cubic Rashba 
Hamiltonian (12311 as 



(26) 



d y = ak x ($kl - fc 2 ), d z 



H k 3 =e(k)+ J2 ^( k ) CT 

i 

where d x = ak y (3k 2 — k 2 ) 
ande(k) = f^-- The eigenvalue associated with the spin 
index fi (fi = ±1) is 



E^k) =e(k) +piak 3 , 



with the eigenfunction 

Vv( r ) = 



fie 



i:\0 



(27) 



(28) 



where As is the area of the system. 

The fc-cubic Rashba model has been widely used to 
study the spin Hall effect in 2DHGr^2&2i however, no 
sufficient attention has been paid to the corresponding 
spin operators. For example, although Hamiltonian ([23]) 
is written in terms of the Pauli matrices a, the a matrix 
is no longer related to the spin directly. The correct spin 
operators in the fc-cubic Rashba model, as described in 
Appendix A, are expressed as 
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Clearly, the coefficient So and the Rashba coefficient a 
have the same dependence on F and L Z1 thus we have 



So = 



am,Q 
h 2 j 2 ' 



(34) 



S z is related to a z , while S x (S y ) consists of two parts: the 
diagonal part linear in k y (k x ) and the non-diagonal part 
quadratic in k±. The diagonal part, which relates the 
wave vector k y (k x ) with S x (S y ), will give the main con- 
tribution to C1SP. The velocity operator in the fc-cubic 
Rashba model can also be obtained by the projection 
technique, 



hk x 3ia 2 2 
v x = h — (fc_(T+ - fc+cr-j, 



(35) 



which is consistent with the relation v x = j-dH/,3 /dk x . 

B. Impurity Vertex correction 

Now, we calculate the ESS in the framework of the lin- 
ear response theory based on fc-cubic Rashba model (|23p . 
In doing this we take the vertex correction of impurities 
into account. The free retarded Green function has the 
form, 



G R (k,E) = 



E - e(k) + djOj 
(E- E+ + ir])(E - + irj) ' 



(36) 



where rj is an infinitesimal positive number. We as- 
sume impurities to be distributed randomly in the form 
V r (r) = V J2i <5( r - R-i)) where V is the strength. With 
the Born approximation, the self-energy, diagonal in the 
spin space, is given by 



Zm[£*(k,£)] 



(37) 



where rii is the impurity density, and the density of states 
for two spin-split branches of the HHI subband reads 



D±(k) 



m h 
2nh 2 



1 ± 



Srrihak 



rr- 



(38) 



So the configuration-averaged Green function is given by 
E - e(k) + iV + J2i dm 



G R (k,E) = 



(E-E, 



ir )(E-E_ 

h 



iTo) 



(39) 



where Tq = — Im{H R (]t, E)] = ^ and r is the momentum 
relaxation time. In the ladder approximation, the Strda- 
Smrcka formula (|f 5| for the ESS \ wn l reduce to 



Xa/3 



n J 2n{ dEjJ (2^ Jr 



(40) 

where S is given by Eqs. (j2"9]) - l|3I[) and the vertex function 
Tp(k) satisfies the self-consistent equation^ 



T/3 = vp + n 4 Vo 2 



d 2 k 
(2nf 



G R (k,E)T f3 G A (k,E). (41) 
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Suppose the electric field is along the x-direction, we solve 
the vertex function iteratively, and get the first-order 



correction to as 
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(42) 



r 



Note that S± and To are independent of 9 and all terms 
in the numerator of the integrand contain something like 
exp(±i9) etc., so the integral over 9 from to 2ir in 
Eq. (|42p vanishes. Furthermore, the higher order terms 
for the vertex correction vanish either, which is quite dif- 
ferent from the vertex correction in the linear-k Rashba 
model.— The same situation occurs for T y . The above 
results agree with the work by Schliemann and Lossi^ 
The calculation of the spin polarization is straightfor- 
ward, and to the lowest order in Fermi momentum k± 
and a, only the term proportional to So contributes to 
the spin polarization. The final result reads 



Xyx 

Xxx 



~Xxy 



eS Tm h E F 

H 3 TT 
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= 0, 



Son h j-, (43) 
(44) 



where rih is the hole density, Ep is the Fermi energy, and 
only the leading term in Ep is retained. 

In the relaxation time approximation the longitudinal 
conductivity of 2DHG equals to 



_ e 2 rEp 

O ' xx — T7f ■ 
/T7T 

Thus, combining the expressions 
the ratio 



(45) 

and l|4"5"l) . we have 
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The formula above can be also obtained from the ex- 
pression of the spin operator (|3H|) and the velocity oper- 
ator (|35[) by neglecting the non-diagonal part in the spin 
operator and the anomalous part in the velocity opera- 
tor, i.e. S y w Sok x and j x ~ ehk x /mh- Obviously, this 
ratio depends only on the material parameters, but not 
on the impurity scattering nor the carrier density in the 
low density limit. Meanwhile, since both the current and 
spin polarization can be measured experimentally, the re- 
lation (|4"6l ) may be invoked to obtain the fc-cubic Rashba 
coefficient a experimentally. 



C. Comparing CISP of 2DHG and 2DEG 

The CISP of 2DHG manifests itselve several features 
different from that of 2DEG. To illustrate this, let's first 



take a look at the CISP of 2DEG. The electric spin sus- 
ceptibility is given by Xyx = 2eTa e m e /H 2 , where m e is 
the effective mass of electron and a e is the linear Rashba 
coefficient. As shown by Inoue et o/4, the vertex cor- 
rection due to the linear Rashba spin splitting is non- 
trivial. With the longitudinal conductivity of 2DEG 
&xx = e 2 rEp / (h 2 it), we find the ratio of spin polarization 
to the current for the 2DEG is 



(4 e) ) 

Ox) 



Xyx 



2Trm e a e 
eE F 



(47) 



Compared with @gj), we find the CISP of 2DEG is in- 
versely proportional to Fermi energy. This means the 
ratio for 2DEG decreases for heavier doping. This dif- 
ferent Fermi-energy dependence stems from the different 
types of spin orientation for 2DEG and 2DHG. 

The spin orientation, which is the expectation value of 
spin operator S for an eigenstate, is given by 

(kfi\S x \kfJ,) — — Sok sin 9 + fik 2 Si sin 9, 
(kp,\S y \kfj,) = S k cos 9 — [ik 2 Si cos 9, 
(kfj,\S z \k(x) = 0, 



(48) 
(49) 
(50) 



for 2DHG, and 



(M4 e) IM 
(Msi e) IM 



—fi smo 
H cos 9, 
0. 



(51) 
(52) 
(53) 



for 2DEG. In the following, we take (S)fc^ as short for the 
spin orientation above. Eqs. (|5Tj) and ((52)) show that spin 
orientation for 2DEG depends on the spin index fx, which 
has opposite values for the two spin-splitting states. But 
for 2DHG, the first term in Eqs. (gSJ) and (gSJ) is inde- 
pendent of the spin index /i. Hence, when k is small, this 
spin-index-independent term will dominate over the k 2 - 
term, leading to the same spin orientation for the hole 
state with opposite fi. This is quite different from the 
electron case. An interesting question may be raised: 
why the holes with opposite \x have the same spin orien- 
tation? In the following, we will analyze this problem and 
try to find the origin of this particular spin orientation 
for 2DHG. 

Let's first have a look at the electron case. Due to 
the spin-orbit coupling and inversion asymmetry, two- 
fold degeneracy of a subband is lifted. For a given k, 
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sin|e-^|i) z 



cos • 



1\ 

2/Z! 



we denote two spin-split states as 
sinf|-i) 2 

where | ± h) z are the eigenstates of a z . It is easy to ver- 
ify that |+) and |— ) have the opposite spin orientation, 
namely (+\ff\+) = 

Similar to 2DEG, two spin-split hole states in the sub- 

i<P \ 3 



± |) can be constructed as |+) = cos |e 
3x and |-) = -sin|e-^||) + cosll - I). By 



space 
sin 1 1 
Eqs. | 

and S y between ||) and 



2/ _ , , = -sinfe-^ID+cosf, a , 
and © , we can verify the matrix elements of S x 



|) vanish, and (ilS^li) = 
(ilS'Ji) = 0. This indicates that in the subspace | ± |), 
any superposition of |±|) will not give rise to the spin ori- 
entation along the x- or y-direction. Thus it is necessary 
to take the higher order perturbation into account, in 
particular the perturbation from coupling between | ± |) 
and | ± |). 

Now we give the outline on the origin of the hole spin 
orientation by the perturbation procedure ( more system- 
atic method can be found in Appendix A). Suppose the 
HH1± states i &hh± can be expanded as 



hh,± 



(i) 

where ± denotes the ith-order perturbed wave func- 
tion. With the basis |n, A^) [Eq. (fT2)) ] and the Oth-order 
term 
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hh,± 
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hh,± 



1,± 



we have the first-order correction as 
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and the second-order correction reads 

,(2) _ lM)(M|K|2,|)(2,l|flt|i 
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|l,-I)(l,-i| J R|2,-|)(2,-f|K|l,-l) 



(-^l,-* _ i)(E 1 i 
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Here E nt x h stands for the eigenenergy of the state 
\n,Xh)- From Eqs. |(HJ) and ([9]), we can see when 
n = 1 the only nonvanishing terms are (1, |) 



and (1,-115^)11 



). Up to the second-order 
perturbation, two types of terms can contribute to 

{^hh,±\Sx(y)\^hh,±}- 

The first type stems from the first-order perturbation 
by the T-operator in the Luttinger Hamiltonian [Eq. (|3J)], 
which couples |1, -±) (|1, §)) to |1, -§) (|1, §)) [the sec- 
ond term in Eq. ({56")) or ((57)) ]. So the matrix clement 

(^hh,+\Sx\^hh,-) equals to 



(0) 
hh,- 



\SJH> 



(1) 

hh,- 



(1) 



(o) 



hh,+ \'~'x \ * hh, 



47T 2 



L 2 Z P .(60) 



It is obvious that the above formula is just the off- 
diagonal element in S x matrix [Eq. f|29|)] with the first 
term in square bracket of Si [Eq. ((55)) ] retained. This 
gives the quadratic-fc dependence of the spin orientation 
shown as the second term in Eq. ((4"8]) . 

The second type comes from joint action of the R in 
in Luttinger Hamiltonian and the asymmetrical potential 
V a [See Eqs. ((55)) and ((59")) ]. The second-order perturba- 
tion contributes to (^h,h,+\Sx\^hh,+) with 



(2) v 
hh,+l 

(1) 



(0) \ 
hh,+l 



-\*hh,+/ 



97r 6 n 2 (37i + 107 2 )(7i- 2 72 )' 



(61) 



This term is just the diagonal element in Eq. (|29|) . which 
leads to the first term in Eq. (148)) and is resposible for 
the identical spin orientation for two spin splitting hole 
states in small k regime. 

The spin splitting between HH± depends on the cou- 
pling between |1,§) and |1, — |) through higher-order 
perturbation. Different from the electron case, the di- 
rect coupling will not cause the x-direction or y-direction 
spin orientation. Instead, it results from the coupling be- 
tween |1,|) (|1,-|)) and |l,i) (|1,-|>). For two LH1 
states, denoted as &ih,±, such coupling will lead to the 
spin orientations of ^ih.+ opposite to ^hh,+i and that of 
^lh - opposite to #hh,-- Thus the total spin orientation 
of the 2DHG is conserved, though ^hh,+ and ^hh- have 
the same spin orientation in the low hole density regime. 



IV. NUMERICAL RESULTS FOR CISP IN 
2DHG 

Based on the calculated eigenstates and eigenenergies 
of the total Hamiltonian ([TJ) , in this Section we will work 
out the spin polarization by using the Bastin version of 
Kubo formula (|14p in the relaxation time approximation. 
Of course the validity of such approximation depends on 
the vanishing vertex correction as mentioned above. 
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Our numerical results with an expanded basis set of N 
basis functions (N is much larger than 8 used in last Sec- 
tion) shows that for a quantum well with infinitely high 
potential barrier, when increasing N, the eigenenergies 
converge to the exact solutions formulated by Huang et. 
al£& very quickly. For example, for the quantum well 
with width L z = 83A, several lowest hole subbands ob- 
tained with N = 20 are almost identical to the exact 
results. Even for N — 8, the dispersion of the lowest 
heavy and light hole subbands is in good agreement with 
the exact results, demonstrating the validity of the trun- 
cation procedure in last Section and Appendix A. Fig. [1] 
plot the dispersion curves and spin splitting of hole sub- 
bands in the quantum well in the presence of an electric 
field. Due to the heavy and light hole mixture effect, 
the energy minimum of the lowest light hole subband, 
marked by B in the Figure, deviates from the T-point 
significantly. 




HH1- 



0.2 0.3 0.4 

k(1/nm) 



FIG. 1: (Color online) Dispersion relation for a quantum well 
with infinite barrier in an electric field. HHH and LH1± 
denote two lowest heavy- and light- hole subbands, respec- 
tively. The parameters for calculation are taken as: the well 
width L z = 83i, the field strength F = 50kV/cm, 71 = 7 
and 72 = 1.9 

For the electric spin susceptibility, we calculate Xyx 
only, because Xxx = Xyy = and Xxy = Xyx as indi- 
cated by Eq. (f2"T|) . After some algebra, we can divide 
ESS in Eq. (|14p into an intra-subband part Xy X an d an 
inter-subband part Xyxi which are expressed respectively 
as 



Xyx 



A.yx 



eh 

eh 
2^ 



AL^(^l4iM(H^|^)^,(62) 

d 2 k 



(27T)* 



$t({kv\Sy\kv')(ki/\v x \kis))A„A u ,. (63) 



the hole subband. In relaxation time approximation, the 
spectral function A v can be expressed as 



A,, 



2tj 



(64) 



where 77 = 



< 20 

> 

S 15 



10 



-10 



■ intra-subband 

■ inter-subband 
both 




0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 

Fermi energy (eV) 



FIG. 2: (Color online) The calculated ESS for the intra- 
subband term (black square line), inter-subband term (red 
circle line), and their sum (green triangle line). The scattering 
induced broadening r\ is taken as 1.65x 10 _5 eF, corresponding 
to the relaxation time t — 2 x lCP n s. The spin polarization 
peak marked with P corresponds to the energy minimum of 
the lowest light hole subband marked as B in Fig. [1] 

A typical curve for the CISP is plotted in Fig. [H The 
main contribution to CISP comes from the intra-subband 
term, which can be understood by Eq. (p3")) . In the 
limit 77 — > 0, the spectral function A v tends to be the 
delta-function 2irS(E — E u ), making the inter-subband 
term A v A v i to vanish except for an accidental degen- 
eracy. Several features in Fig. [5] are worth pointing out. 
First, in low doping regime where only HH1± states near 
r point are occupied, spin polarization exhibits a linear 
dependence on the Fermi energy. Second, with the hole 
density increased, the spin polarization increases at first, 
then decrease after reaching a maximum value, and even 
changes its sign when the hole density is large enough. 
Third, when the doping is so heavy that the light hole 
subband is occupied, a sharp peak for the spin polariza- 
tion may be observed as marked as P in Fig. O 

To understand these features, we turn back to Eq. (|62|) . 
as main contribution to the spin polarization stems from 
this intra-subband term. Based on numerical results as 
well as Eq. (|4*g|) , we adopt a function J v (k) to express 
the amplitude of the spin orientation associated with the 
subband v. i.e. 



Here 5i denotes the real part, and v and v' stand for 



(Sy 



J v {k) cos( 
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Then, with 



. . ldE v ldE u {k) 



and 



we rewrite Eq. 



,o 47TT „. 

A 2 U = —S(E F - E v ), 



as 



Xyx — . ^ J]] k„J u (k^), 



(65) 



where fc,, is the Fermi momentum with the hole subband 




0.2 0.4 0.6 

k (1/nm) 



1.0 



In the fc-cubic Rashba model, in which only the lowest 
heavy hole subband HH1± is concerned, up to the first- 
order in a, the Fermi momentum can be expressed as 
k F = V2mHE F _ n 2am l E F . Combined with Eq.®, we 
obtain 



Xyx 



nh 3 



SermiaSiEp 

5rF ' 



(66) 



< 

* 

> 



LU 




0.045 0.048 0.051 

Fermi energy (eV) 




0.4 0.6 

k (1/nm) 



FIG. 4: (Color online) The magnitude of the spin orientation 
for the lowest heavy hole subband (a) and the lowest light hole 
subband (b). In (a), the black solid line and red dashed line 
represent Jh+ and Jh-, respectively; while the green dotted 
line and blue dashed dotted line denote the spin- indepen- 
dent part J h and dependent part J h , respectively. The same 
notions are also applied to (b). 



FIG. 3: (Color online) Numerically calculated ESS as func- 
tions of Fermi energy (black-square-line) compared with the 
analytical results (red-circle-line). 

The first term on the right hand side of Eq. (f6"6"]) . re- 
sulting from the spin-independent part, is identical to 
Eq. (|4"5)) ; while the second term, proportional to E\ , can 
be safely ignored in the low density regime. As shown 
in Fig. [31 the analytical results of the electric spin sus- 
ceptibility (Eq. I66[) agree well with the numerical ones, 
demonstrating the applicability of /c-cubic Rashba model 
(|23p in low doping regime. However, for higher hole den- 
sity, numerical results show a drop of the \ due to the 
heavy and light hole mixing effect, which is certainly be- 
yond the simple fc-cubic Rashba model. 



For numerical results, similar to the derivation above, 
we may divide J„ into a spin-dependent part and a spin- 
independent one, namely, J„ M = Jl + fiJ^. Then the ESS 
can be expressed as 

Xyx = X\ x +Xy X , (67) 

in which the spin- independent and dependent part re- 
spectively reads 

*vx - ^rX 3 « — 2 — ' (68) 

4. = jaE**^- («> 

Obviously. y£ depends on the average of Fermi 
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wavenumbers, while Xy X depends on the Fermi wavenum- 
ber difference between two spin-split branches. In most 
cases, owing to the fact that the spin splitting is small 
compared with the Fermi energy, \\ x w ^ dominate the 
spin polarization. In Fig. [4ja), we plot the magnitude 
of spin orientation associated with the subband HH1±, 
denoted by Jh±, and the corresponding spin- dependent 
part, J^, and independent part J l h . They are related 
through J* = [J h+ - J h _) /2 and J l h = {J h+ + J h J) /2. 
Fig. 2] indicates that for most values of k is larger 
than J l h . Compared to the intra-subband contribution in 
Fig. [21 the spin-independent magnitude of the spin po- 
larization J\ [green-dotted line in Fig. Ufa)] has similar 
behavior: first increasing linearly with k, then decreasing 
with k increased, and even changing the sign for larger 
k. 



0.080 



> 
LU 



E r 








LH1 + / 


/.H1- 


E° 
























k" 


k F1 

1 


4 


!k° 


C 


k F2 



0.1 0.2 

k(1/nm) 



polarization cancel each other to some extent, resulting 
in 

Xv* = ^P?- + k? + )Jf + (kf + - kl)J?}. (71) 



kl 



As J\ is much smaller than jf or Ji±, and kf, 
both terms in Eq. l|7"Tj) are small compared to the case 
when only LH_ is occupied. Apparently, the peak width 
depends on the spin splitting between LH- and LH+. 

The temperature dependence of the peak is plotted in 
Fig. [6l Near the polarization peak, if we only take into 
account LH1±, ESS is expressed by 



(72) 



At zero temperature, the Fermi distribution function 
f(E) becomes the step-function 0(Et — E), which re- 
produces the above analysis. At finite temperature T, if 
we approximate kf^Ji^k®^) ~ ^kfJi, and expand the 
Fermi distribution function at large ksT as f(E) — 
|(1 — ) ( is Boltzmann constant), then Eq. ([72"l) 

reduces to 



Xyx 



erkfJi El - El 
2nh -ik B T 



(73) 



So ESS is proportional to the ratio of the spin splitting 
of the LH1 subband, Ef + — Ef_, to thermal energy fc^T. 
When ksT is much larger than the spin splitting, this 
pronounced spin polarization peak will smear out. 



FIG. 5: (Color online) Dispersion relation of the lowest light 
hole subband LH1±. 

A pronounced peak of CISP may appear when the 
Fermi energy just crosses the bottom of the lowest light 
hole subband LH1-. As amplified in Fig. O in the dis- 
persion relation of the subband LH1±, the wave num- 
bers corresponding the energy minimum Ef ± deviate 
from the k — point significantly. Around the energy 
minimum the energy dispersion can be approximated as 

E llM (k) = E^ + § a2 y } (fc ~ Assuming the above 

energy dispersion and a constant magnitude of Ji^, we 
obtain 



< 


20 -, 
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(hbar/ 


15- 
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spin 
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lectr 
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LU 






Xyx — 2~ T^lfl,Jlfl(^lfJ,)i 



(70) 



0.076 0.078 0.080 0.082 
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FIG. 6: (Color online) The spin polarization peak at three 
different temperatures, 1.2K, 12K and 120K. 



where k° ~ (k^ 1 +k^ 2 )/2, and fc^ 1 and k± 2 respectively 
denote two different Fermi wave numbers for LHl± (Fig. 
[5]). By Eq. (ITO)) and Fig. H[b), we can see since J; + and 
Ji^ are large in the absolute value but almost opposite in 
the sign, when Ef, > Ep > Ef_, a large spin polarization 
eTkf_Ji-/(2irh) is expected; on the other hand, when 
Ep > Ei + > E®_, the contributions of LH1± to the spin 



Now let's estimate the magnitude of the averaged 
CISP. In the fc-cubic Rashba model with an applied filed 
F = 50kV/ cm, Eq.® gives S = 2.74A for L z = 83 A 
and So = 5.77A for L z = 100A If typical relaxation 
time t is taken to be 2 x lCP n s and an in-plane electric 
field strength Eq = lOV/cm, the Fermi sphere will be 
shifted by Ak = eE a r/h = 3x lCT 3 ! -1 . Substituting 
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020- 




0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 



Fermi energy (eV) 



energy, numerical results show that the spin polarization 
may be suppressed and even changes its sign. We predict 
and explain a pronounced spin polarization peak when 
the Fermi energy crosses over the subband bottom of the 
LH-. We also discuss the possibility of measuring this 
spin polarization peak. 
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FIG. 7: as functions of the Fermi energy. The inset: 

(s ) 

1 - !LL as functions of applied field strength F. 



APPENDIX A: DERIVATION OF THE fc-CUBIC 
RASHBA HAMILTONIAN 



the above data into Eq. (|43"|) , we obtain (S y ) /rih = 0.831% 
for L z = 83 1, and (S y )/n h = 1.75% for L z = lOOi. 
Since So is proportional to L\, the spin polarization is 
very sensitive to the thickness of quantum well. Hence, 
it is preferable to detect the CISP in a thicker quantum 
well experimentally. The above estimation gives the same 
order of magnitude for the spin polarization observed in 
Silov's experiment^. In Fig. [Jj we plot the averaged 
spin polarization (S y ) jnh as functions of the Fermi en- 
ergy and functions of the field F in the inset. The CISP 
is saturated about 2% when the field is enhanced. 



SUMMARY 



In this Appendix, we present the detailed derivation 
of the fc-cubic Rashba model by means of the per- 
turbation method.-^^L^^Mi^ First we truncate 
the Hilbert space of the basis wave functions (fl"2"]) into 
the subspace with only the lowest eight states Go = 
{\n,Xh),n = 1,2; Xh = ±|,±|}. As described in the 
Sec. II, by comparing the lowest HH and LH subband 
dispersion with the exact solution, the accuracy of such 
truncation procedure has been verified. The truncated 
subspace Qq can be further cast into two sub-groups, 
Qi and Qi. Qi contains two lowest heavy hole states 
{| 1,3/2), |1, — 3/2)}, while Qi keeps the other six states, 
{|1, 1/2), |1, -1/2), |2,3/2), |2,-3/2), |2, 1/2), |2, -1/2)}. 
In this case, the Hamiltonian in the subspace Qq can be 
written in the form of block matrices as 



In conclusion, we have systematically investigate the 
current induced spin polarization of 2DHG in the frame 
of the linear response theory. We introduce the physical 
quantity of the electric spin susceptibility \ to describe where 
CISP and give its analytical expression in the simplified 
/c-cubic Rashba model. Different from the 2DEG, the 
CISP of 2DHG depends linearly on the Fermi energy. 
The difference of CISP between 2DHG and 2DEG re- 
sults from the different spin orientations in the subband 
of carriers. We propose that /c-cubic Rashba coefficient 
of 2DHG can be deduced from the ratio of spin polariza- 
tion to the current, which is independent of the impurities 
or disorder effect up to the lowest order. We have also 
carried out numerical calculations for the CISP. The nu- 
merical results are consistent with the analytical one in 
low doping regime, which demonstrates the applicability 
of fc-cubic Rashba model. With the increase of Fermi and 



#8x8 — 



#? 



#6x2 - # 2x6 



#2x2 #2x6 
#6 x 2 #6 x 6 



P(l) 
P(l) 



T \ 

Tt 

eFG(2,T) 

eFG(2,l) 

R(2,l)k+ 

-P(2,l)fc_/ 



(Al) 



(A2) 



(A3) 
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/ Q(i) 





P(l,2)fc+ 





eFG(l,2) 





\ 










-P(l,2)fc_ 





eFG(l,2) 


R(2,l)k- 





P(2) 








T 







-P(2,l)fc+ 





P(2) 


T t 







eFG(2,l) 








T 


0(2) 







V o 


eFG(2,l) 


rt 








0(2) 


/ 



(A4) 



Here P(n), Q{n), G(n, m), m) are given by 
ft 2 



P(rc) 
Q(n) 



2mo 
2mo 



(7i+72)fc 2 + (7i- 272)(^) 2 
(7i-72)A: 2 + (7i + 272)( 1 -) 2 



(A5) 
(A6) 



4£ z nm((-l)"+ m - 1) 
G(n,m) = 2T2 ' ( A7 ) 



^ft 2 73 2mm((-l)"+ m - 1) , 4 N 
R(n, m) = -2V3—^ (A8) 

Our aim is to perform a transformation which decou- 
ples the groups Gi from Q 2 , i.e. to make the off-diagonal 
part #2x6 and #6x2 vanish up to the first-order in k and 
F. We divide the total Hamiltonian (|A1[) into three parts 



H 



8x8 



# + #i+# 2 . 



(A9) 



The first term #0 is the diagonal matrix elements of 
# 8 x8, given by 



#0 



# 



(o) 

6x6 



(A10) 



with H^ 2 = Diag[P(l),P(l)} and H^% 
Diag[Q(l), Q(l), P(2), P(2), Q(2), Q(2)]. 
The second term #1 is given by 



Hi 







# 



(i) 

6x6 



(All) 



where #6x6 = H ex6 — H^ 6 . The third term # 2 contains 
the non-diagonal part H 2x6 an d #6x2 



#2 = 



#2x6 
#6x2 



(A12) 



There are three types of perturbation terms in #1 and 
#2: (l)The k-linear R term couples the state \n, §) 
(|n, — I)) with |m, 5) (|m, — |)), where n and to must 
be of opposite parities due to the presence of k z = —id z ; 
(2) The k-quadratic T term couples \n, §) (\n, — |)) with 
|n, — i) (|n, ^)); (3) The asymmetric potential V a couples 
the states with the same spin index and different parities. 



The perturbation procedure is as follows. First H 2 will 
be eliminated by the canonical transformation as 

r(l) 



#8x8 = cxp[-(7 (1) ]# 8X 8exp[(7 (1) ] 

= # 8 x8 + [#8X8, U (1) ] + ^[[#8X8, U^},U ( ^} 



(A13) 



in which [A 1 ) is chosen such that 

# 2 + [Hq, # (1) ] = 0, 
and the matrix elements read 



#. 



(1) _ (H 2 )af3 

al3 ~ E a — Ep" 



(3, 



(A14) 



where E a denotes the energy of the band a at the T 
point (k=0). After the canonical transformation, the new 
Hamiltonian is given by 

#« 8 = #„ + #! + \[H 2l #«] + [#!,#«] 

+l[[H 1 ,U^] ) U^} + ---. (A15) 

The # , #1, i[# 2 ,C/ (1) ] and i[[#i, #W], #«] have 
the block-diagonal form, while [Hi,U^] is non- block- 
diagonal and contains new terms first-order in k. So we 
divide #gx8 i n t° three parts again 



#. 



(i) 

8x8 



#0 + # 1 (1) +#W, 



(A16) 



in which #j X) = #! + i[# 2 ,£/ (1) ] + i[[#r, Ul%UW], 

and #2 1 ' ) = [#1, U^]. We perform the second canonical 
transformation U^- 2 \ given by 



(2) 

a/3 



(H 2 ^) a 

E a — Ep' 



(3. 



(A17) 



This makes the non-diagonal block matrix #2 zero, 
leading to the Hamiltonian 



#8X8 - #0 + # 



[ 1) + 1 -[HP,U^] + [Hi 1 \u^] 

+i[[# 1 1) ,t/«],[/«] + ---.(A18) 
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(2) 

Now the non-block-diagonal terms of -ffg x8 vanish up to 
the desired order in k and F. Finally, by mapping the 

(2) 

Hamiltonian -ffg x8 into the lowest heavy hole subbands, 



we obtain the £;-cubic Rashba Hamiltonian Eq. (|23|) . 

To obtain the corresponding spin operators in the low- 
est heavy hole basis, we should apply the same canon- 
ical transformations and to the spin opera- 
tors Si (i = x,y,z). In the 8-state subspace Go, we 
find that the spin operator has the block-diagonal form 

Si = Diag[Si ,S^] (* = x tUi z )i because there are no 
matrix elements between the states with different con- 



finement quantum number n. Therefore S, 
matrix, given respectively by 



S~. 
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is a 4 x 4 



(A19) 



(A20) 



where A 



„SvSv } 
'41 



„8vSv 
'41 



is a parameter as already given 



by Winkler for several materials^ and F is the field 
strength. If we neglect other asymmetrical potentials and 
only consider the Rashba term, then the total Hamilto- 
nian is H = Hl + V c + Hr. Applying the same perturba- 
tion procedure as in the appendix A, we find that both 
the Hamiltonian and the spin operator have the identical 
structure to the asymmetrical potential case, as well as 
the same effective mass run, S\ and expression Eq. (j34"| . 
except for the Rashba coefficient given by 



3AL 2 
4tt 2 ' 



and the spin operator parameter 



(B2) 



So = 



3Xm L 2 z 



(B3) 
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Then we apply the transformations and to 
the spin operators, obtaining the new spin operators Si = 
Si + [Si, [/W] + [Si, £/ (2) ] as presented in Eqs. @, §U§ 
and JED). 



APPENDIX B: HOLE RASHBA TERM 

The hole Rashba term has recently attracted many 
researchers' attentions. 32 ' 36 i 47 The hole Rashba term 
breaks the inversion symmetry, 35 ' 37 and is expressed as 






2 K ~ 





2 K + 





ik- 
















iV3h. 
2 ft 






iV3 



(A21) 

The hole Rashba coefficient a here is proportional to L z , 
while for the asymmetrical potential case it depends on 
L z . So in most realistic quantum wells, the contribution 
from the asymmetrical potential plays more important 
role than the hole Rashba term, at least one or two or- 
ders of magnitude larger. The physical reason for this 
may be understood from the origin of the hole Rashba 
term. The more general form of the Hamiltonian should 
be H — H^ p + V c + eFz, where the multi-band k • p 
Hamiltonian H^.p includes not only the heavy and light 
hole bands, but also the conduction band, spin split-off 
band and remote bands. When we project the Hamilto- 
nian into the subspace of the heavy and light hole bands, 
the combined effects of the eFz and k ■ p mediated by 
other bands lead to the hole Rashba term, which has 
much smaller influence than that coupled by the asym- 
metrical potential directly. Therefore, hole Rashba term 
, (Bl) is neglected in the present article for simplicity. 
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